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Abstract 

We describe the automatic synthesis of a global nonlinear con¬ 
troller for stabilizing a magnetic levitation system—a simplified model 
for the German Transrapid system. A systematic state-space design 
method for determining the global switching points of the controller 
is presented. The synthesized control system can stabilize the maglev 
vehicle with large initial displacements from an equilibrium, and pos¬ 
sesses a much larger operating region than the classical linear feedback 
design for the same system. 

The controller is automatically synthesized by a suite of compu¬ 
tational tools that visualize and model the state-space geometry and 
topology of a given system, use a novel technique of “flow pipes” to 
plan global reference trajectories in state space, and navigate the sys¬ 
tem along the planned trajectories. This work demonstrates that the 
difficult control synthesis task can be automated, using programs that 
actively exploit knowledge of nonlinear dynamics and state space and 
combine powerful numerical and symbolic computations with spatial- 
reasoning techniques. 

Keywords. Magnetic levitation; control system synthesis; state- 
space methods; artificial intelligence; numeric/symbolic processing; 
computer-aided design. 
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1 Introduction 


Magnetically levitated trains provide a high speed, very low friction alterna¬ 
tive to conventional trains with steel wheels on steel rails. Several experimen¬ 
tal maglev systems in Germany and Japan have demonstrated that this mode 
of transportation can profitably compete with air travel. More importantly, 
maglev transportation can ease traffic congestion and save energy [2, 3, 5]. 

The maglev transportation uses magnetic levitation and electromagnetic 
propulsion to provide contactless vehicle movement. There are two basic 
types of magnetic levitation: electromagnetic suspension (EMS) and electro¬ 
dynamic suspension (EDS). In EMS, the guideway attracts the electromag¬ 
nets of the vehicle that wraps around the guideway. The attracting force 
suspends the vehicle about one centimeter above the guideway. In contrast, 
the EDS system uses repulsive force, induced by the magnets on the vehicle, 
to lift the vehicle. 


2 The Maglev Model 


An attractive system such as the EMS system is inherently unstable. We 
consider the control design for stabilizing an EMS-mode train traveling on 
a guideway—a simplified model for the German Transrapid experimental 
system. The Transrapid system is schematically shown in Figure 1. It uses 
attractive magnetic forces to counterbalance gravitational forces. 

The state equations for the magnetically levitated vehicle and the guide¬ 
way are described by 


( dx z(Vi — Rx) xy 

dt LqZq z 

dy LqZqx 2 


~t= 9 


2 mz 2 


n = y 


( 1 ) 


where the state variables x, y, and z represent coil current, vertical veloc¬ 
ity of the vehicle, and vertical gap between the guideway and the vehicle, 
respectively. The control parameter is the coil input voltage Vi. The other 
parameters are the mass of the vehicle m, the coil resistance R, the coil in¬ 
ductance Lo and the vertical gap Zq at the equilibrium, and the gravitational 
acceleration g. Details of the derivation of the model are discussed in [4]. 
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Electromagnetic Maglev 




LIFT 


(b) 


Figure 1: EMS maglev system for high-speed ground transportation, repre¬ 
senting a simplified drawing of the German Transrapid design, (a) Electro¬ 
magnetic suspension; (b) Detail of a suspension magnet, superimposed on 
the field distribution. 
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The system has one equilibrium state at which the magnetic force exactly 
counterbalances the force due to gravity and the vehicle has no vertical veloc¬ 
ity and acceleration. However, the equilibrium is a saddle node which is not 
stable. The control objective, therefore, is to stabilize the vehicle traveling 
down the guideway and maintain a constant distance between the vehicle and 
the guideway despite any roughness in the guideway. The available control 
input is the coil input voltage Vi, as in the model (1). We further assume 
that Vi is produced by a buck converter capable of delivering any voltage 
from 0 to 300 volts. 

A linear control design for the maglev system described in [4] uses the 
pole-placement method. The system is first linearized around the equilib¬ 
rium. The linearized system has unstable poles, i.e., the poles in the right-half 
of s-plane. A linear feedback is introduced to move the poles to the desired 
locations in the left-half of the s-plane. Such a control design can bring the 
system back to the equilibrium with an initial displacement of 0.2mm from 
the equilibrium. The linear controller saturates at the beginning for larger 
initial displacements. This is because the linearized model no longer approx¬ 
imates the original system well in regions far away from the equilibrium. A 
nonlinear control law such as a bang-bang control that respects the nonlin¬ 
earities of the system is, therefore, necessary to proceed the linear feedback 
control. However, the real challenge for the nonlinear design is to determine 
the global control switching points. 

3 State-Space Control Trajectory Design 

We demonstrate a nonlinear control design—a switching-mode control—for 
the maglev system with large initial displacements from the equilibrium. 
We also present an automatic method for designing such a controller. The 
nonlinear controller brings the system to the vicinity of the equilibrium and 
then switches to the linear controller. 

For the purpose of demonstration, we assume that the vehicle is displaced 
from the equilibrium in the direction further away from the guideway. We 
will concentrate on the global design of the control reference trajectories and 
assume that a linear feedback design is available as soon as the system enters 
the capture region of the linear controller. 
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3.1 Modeling state-space geometry and topology 

The global control law is designed by exploiting the nonlinear dynamics in 
state space. We have developed a suite of computational tools for the de¬ 
sign and analysis of nonlinear control systems. The two main programs, 
MAPS and Phase Space Navigator, work together to visualize and model 
the state-space geometry and topology of a given system, plan global refer¬ 
ence trajectories in state space, and navigate the system along the planned 
trajectories. MAPS computationally explores the state space of a dynamical 
system and characterizes the state space with stability regions and trajectory 
flow pipes that bundle together trajectories with the same qualitative fea¬ 
tures [6]. Phase Space Navigator uses the state-space characterizations for 
different control paxameter values to synthesize a composite reference tra¬ 
jectory, consisting of path segments connected at intermediate points where 
control parameters change [7]. 

We ran MAPS on the maglev model. The state variables x, y, and z in 
the model are scaled by 1, 10 3 , and 2 x 10 4 , respectively. The parameters of 
model are assumed to be: L 0 = 0.1 h, z 0 = 0.01m, R = lfl, m = 10000% 
and g = 9.8m/sec 2 , typical of a large vehicle lift magnet. Assume at the 
equilibrium, the power supply delivers 140 volts, i.e., VJ = 140. MAPS 
explores in a region of the state space bounded by the box {(#, y, z)\x € 
[0,400], y G [—300,350], z G [0,600]} and finds the following equilibrium 
point: 

saddle: #(140. 0. 200.) 

eigenvalues: -17.004+22.963i 
-17.004-22.963i 
24.007 

eigenvectors: #(.23604 .97174 0) 

#(.51331 -.55588 .65384) 

#(.30157 .73255 .61027) 

With the information about the stable eigenvectors of the saddle, MAPS 
computationally determines the stable manifold (often called separatrix) of 
the saddle. MAPS uses a set of trajectories evenly populated on the stable 
manifold to approximate the manifold. The set of trajectories are obtained 
by backward integrations starting at initial points in the plane spanned by 
the stable eigenvectors within a small neighborhood of the saddle. The stable 
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Figure 2: Stable and unstable manifolds of the saddle for Vi = 140 (yz- 
projection). 

manifold is geometrically modeled by a tessellation on these trajectory points 
with a triangulation method. 

Figure 2 shows the trajectories on the stable manifold and along the 
unstable eigenvectors of the saddle in the yz-projection of the state space. 
The stable manifold of the saddle separates the state space into two halves: 
trajectories in the upper-half approach z —* oo along one of the unstable 
eigenvectors, corresponding to the case in which the vehicle falls off the rail; 
and trajectories in the lower-half approach z — 0 plane along the other 
unstable eigenvector, corresponding to the case in which the train collides 
with the rail. 


3.2 Synthesizing a global stabilization law 

For an initial displacement above or below the equilibrium, the uncontrolled 
system will follow either a trajectory traveling upwards with increasing z and 
leaving the bounding box or a trajectory traveling downwards and hitting the 
z = 0 plane. To stabilize the system at the equilibrium, it is necessary to 
synthesize a new vector field on both sides of the stable manifold so that 
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Figure 3: Stable and unstable manifolds of the saddle for Vi = 300 (yz- 
projection). 

trajectories travel towards the stable manifold of the saddle. We consider 
only the top-half here. 

Phase Space Navigator first explores the state space of the model for 
different values of Vi and discovers that the larger the Vi is, the further away 
the stable manifold is from the z = 0 plane. For V{ = v > 140, the region 
sandwiched by the stable manifold of V = 140 and that of Vi = v has the 
desired property—the vector field of Vi = v in this region is pointed towards 
the stable manifold of Vi = 140. When v — 300, the region is maximized. 

Similarly, Phase Space Navigator finds that the model with = 300 has 
a saddle node at (300.,0., 428.57). The stable manifold of the saddle has a 
similar structure as that of Vi = 140 case, but with larger z coordinate. The 
yz-projection of the stable manifold and unstable trajectories for V = 300 is 
shown in Figure 3. 

Figure 4 shows the yz-projection of the sandwiched region discussed 
above. The region is bounded by three pieces of triangulated surfaces, the 
top boundary shown in Figure 5(a), the side one in Figure 5(b), and the 
bottom one in Figure 5(c), respectively. The top boundary represents the 
stable manifold of the saddle for Vi = 300 and the bottom one for Vi = 140. 
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The trajectories of Vi = 300 flow into the region from the side boundary 
in Figure 5(b) and leave the region at the bottom boundary in Figure 5(c). 
There is no flow across the top boundary in Figure 5(a). 

The Phase Space Navigator determines that the above region consists 
of two collections of trajectories: the bundle with Vi = 300 that flow into 
the region on the side boundary and pierce through the bottom one, and 
the bundle with Vi = 140 on the bottom boundary that approach the de¬ 
sired equilibrium asymptotically. Each bundle of the trajectories is called 
a flow pipe that groups together a collection of trajectories exhibiting the 
same qualitative behavior. (For more details about the flow pipes, see Ref¬ 
erence [7].) With the flow pipe characterization of the state-space trajectory 
flows, the Phase Space Navigator searches for control trajectories in this set 
of flow pipes and finds a sequence of flow pipes that lead to the desired 
goal: the composite of the trajectory bundles with Vi = 300 and Vi = 140, 
glued together at the stable manifold of Vi = 140 represented by the bottom 
boundary. As a result, all the trajectories of Vi = 300 within the region 
can be brought to the equilibrium by switching to Vi = 140 as soon as the 
trajectories hit the bottom boundary. 

Since the bottom boundary is an approximation to the true stable mani¬ 
fold, the trajectories of Vi = 140 on the boundary can only get close to the 
desired equilibrium. The closeness depends on the quality of the manifold 
approximation. As the trajectories enter a small neighborhood of the equi¬ 
librium, we use a linear feedback controller, such as the one discussed in [4], 
to stabilize the system at the equilibrium. Figure 6 and 7 show the syn¬ 
thesized control reference trajectories originating at different initial displace¬ 
ments from the equilibrium: 1mm, 4mm, 4.5mm, and 5mm. The controller 
is able to recover from the first three initial points that are within the region. 
The last point is outside the region and thus uncontrollable. 

3.3 Evaluating the control design 

The synthesized global control law is a switching-mode one that changes 
parameters at the switching surface—the stable manifold of Vi — 140. It is 
able to bring trajectories originating from any states within the controllable 
region to a local neighborhood of the saddle. 

The responses of the controller with respect to the different initial dis¬ 
placements are shown in Figures 8. For all the controllable initial displace- 


8 




Figure 4: The controllable region in yz-projection. 
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(c) 

Figure 5: The boundaries of the controllable region in t/z-projection: (a) top 
boundary; (b) side boundary; (c) bottom boundary. 
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Figure 6: The synthesized control reference trajectories originating from four 
different initial states, together with the controllable region in yz-projection. 
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Figure 7: The synthesized control reference trajectories originating from four 
different initial states, together with the controllable region in zx-projection. 
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ments, the controller is able to bring the system as close to the equilibrium 
as within less than 0.2mm in displacement—a distance well captured by the 
linear feedback controller. 

By exploring the state-space geometries of the maglev system, our pro¬ 
grams are able to automatically determine the switching points for the global 
controller. The global controller has significantly enlarged the operating re¬ 
gion of the linear feedback controller. With the geometric representation of 
the controllable region in state space, the programs precisely determined that 
the maximum recoverable displacement is 4.55mm. Our simulation geomet¬ 
rically explained the observation in [4] that the vehicle would fall off the rails 
with a displacement of 5mm or larger. 

Many issues remain to be addressed in order to make our control design 
practical. Future work includes the optimization of the design with respect 
to response time and the study of the effect of uncertainties in the model and 
noises in the environment on the design. 

4 Conclusion 

We have designed a high-quality global nonlinear controller for stabilizing 
an EMS-mode maglev vehicle traveling above guideway. The synthesized 
control reference trajectory consists of a sequence of trajectory segments, 
connected at intermediate points where control voltage switches. At run¬ 
time, the controller will cause the system to track the reference trajectory 
and reactively corrects local deviation from the desired trajectory. We have 
illustrated a systematic method for designing the global switching points of 
the nonlinear controller. The simulation showed that our design has a much 
larger maximum recoverable displacement than the classical linear feedback 
design. This work demonstrates that the difficult control synthesis task can 
be automated, using programs that actively exploit knowledge of nonlinear 
dynamics and state space and combine powerful numerical and symbolic 
computations with spatial-reasoning techniques. 
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Appendix A 


The equilibria for V,- = 140 and V,- = 300: 


; Vi = 140 
((#(140. 0. 200.) 

(-17.004+22.963i 
-17.004-22.963i 
24.007) 

(#(.23604 .97174 0) 
#(.51331 -.55588 .65384) 
#(.30157 .73255 .61027)) 
saddle)) 

; Vi * 300 
((#(300. 0. 428.57) 

(-21.411+21.394i 
-21.411-21.394i 
21.394) 

(#(.54819 .83636 0) 
#(.43085 -.65947 .616) 
#(.23229 .71052 .66423)) 
saddle)) 


equilibrium point #(x y z) 
eigenvalue 1 
eigenvalue 2 
eigenvalue 3 
eigenvector 1 
eigenvector 2 
eigenvector 3 
stability type 

equilibrium point #(x y z) 
eigenvalue 1 
eigenvalue 2 
eigenvalue 3 
eigenvector 1 
eigenvector 2 
eigenvector 3 
stability type 
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Appendix B 


This appendix contains the initial, switching, and final points of the syn¬ 
thesized reference trajectories starting from four different initial states, as 
shown in Figure 6 and 7. Each point in the print-out is a vector of ( t , x, y, z). 

; trajectory 1 

(init-pt #(0. 140. 0. 220.)) 

(switch-pt #(.013441 162.98 4.4416 221.42)) 

(last-pt #(.12704 139.39 .078872 201.75)) 

; trajectory 2 

(init-pt #(0. 140. 0. 280.)) 

(switch-pt #(.10156 276. -74.214 338.79)) 

(last-pt #(.32041 139.26 -.19212 196.51)) 

; trajectory 3 

(init-pt #(0. 140. 0. 290.)) 

(switch-pt #(.19536 264.69 -105.03 340.68)) 

(last-pt #(.42865 142.4 3.2661 202.69)) 

; trajectory 4 

(init-pt #(0. 140. 0. 300.)) 

(init #(0. 140. 0. 300.) is outside controllable region) 
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